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Several kinds of measurements are combined in an attempt to obtain a consistent estimate 

o ■ 

■ of the spectrum and composition of the primary cosmic radiation through the knee region. 

Assuming that the knee is a signal of the high-energy end of a galactic cosmic-ray population, 
, I discuss possible signatures of a transition to an extra-galactic population and how they 

I/-) ■ might be detected. 

<n : 

> ■ I. INTRODUCTION 

m : 
in . 

The cosmic-ray spectrum extends from the sub-GeV region to at least 10 11 GeV (Fig.ffl) 
Up to 100 GeV and somewhat higher, measurements with magnetic spectrometers flown 
above most of the atmosphere provide good momentum resolution along with identifi- 
cation of the charge and mass of individual primaries. Measurements with calorimeters 
continue to identify individual primaries to beyond 100 TeV, but with larger systematic 

& ! 

uncertainties in the energy assignment. In the PeV region and beyond, the cosmic-ray 
| intensity is too low for direct measurements; only indirect measurements of air showers 

•rH . 

^ ■ from the ground are possible. Since the particles are not identified on an event-by-event 

basis, the energy spectrum derived from measurements of air showers is given as an "all- 
particle" spectrum, in terms of energy per particle rather than energy per nucleon. In the 
air-shower regime, identification of the primary mass is made in one of several indirect 
ways on a statistical basis, complicating the search for features in spectra of individual 
elements. 

The gyroradius of a proton in a typical galactic magnetic field of 3 ^Gauss is about 
half a parsec at one PeV. The parsec scale is also typical of the size of structures in the 
interstellar medium driven by supernova explosions. For particles with higher energy and 
larger gyroradii, diffusion in the interstellar medium may become less efficient. Estimates 



Energies and rates of the cosmic-ray particles 
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FIG. 1: Spectrum of cosmic rays 

of the maximum energy of particles accelerated at supernova shocks are around the same 
energy. It is therefore natural to guess that the knee of the cosmic-ray spectrum around 3 
PeV has something to do with the end (or at least the beginning of the end) of the galactic 
cosmic-ray spectrum. 

Peters [l| described the consequences for energy dependence of the primary composition 
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if the spectrum is characterized by a maximum rigidity, R c , which could be associated 
either with propagation or with acceleration (or both). The relation between rigidity and 
total energy is 

Pc 

R = (i) 

Z e 

where Fc ~ E tot * s the- total energy of a nucleus of charge Z and mass A. If cosmic 
rays are classified by energy per particle, as is the case for air shower measurements, then 

the spectrum should steepen first for protons, then for helium, then for the the CNO 

[ 

group etc. Elsewhere [2| I have called this sequence the "Peters cycle". Several air shower 
measurements show some evidence that the spectrum becomes progressively enriched in 
heavy nuclei through the knee region. The clearest evidence for the Peters sequence for 
several groups of nuclei comes from analysis of the KASCADE experiment |3j. 



II. COMPARISON OF DIRECT AND INDIRECT MEASUREMENTS 

Data from direct measurements above 100 TeV are sparse, leaving a gap that is bridged 
by emulsion chamber data with low statistics below a PeV and by the threshold region of 
air shower experiments above 100 TeV. Apart from the technical problems of low statistics 
and systematic threshold effects, there is also the problem that air shower experiments do 
not identify individual primary nuclei. In addition, the efficiency of ground arrays depends 
strongly on primary mass in the threshold region. Nevertheless, it is possible to form a 
fairly consistent picture of the cosmic-ray spectrum up to the knee. 

Figure El shows several measurements of the primary spectrum through the knee region. 
Up to 100 GeV there are good measurements of the spectra of protons and helium with 



magnetic spectrometers flown in spacecraft |4J and high-altitude balloons 



Using the 



P| and high-altitude balloons P, (| . 
AMS 141 and BESS 5fl measurements, together with measurements of heavier nuclei at 

n 

10 GeV/nucleon .7], I have converted the spectra of protons, helium, CNO, Ne-Si and Fe 
from energy per nucleon to total energy per particle and combined them to give the all 
particle spectrum at low energy (shown as the solid line in Fig. 0). 
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FIG. 2: All-particle spectrum. Upper panel illustrates an extrapolation from direct measurements 
at low energy to 1 PeV. Lower panel shows detail of air shower measurements in the knee region. 
The two lower sets of points are offset respectively by 1/10 and by 1/100 for clarity (see text). 

The highest energy direct measurements in which individual nuclei are identified di- 
rectly are from emulsion chamber experiments. Measurements of protons and helium from 
n ri 

RUNJOB |8| and JACEE |9| are shown in Fig. EH- Preliminary data on protons and he- 
lium from the ATIC thin ionization chamber [llj (not shown here) are consistent with 
a smooth power-law extrapolaton between the BESS and AMS spectromter data at low 
energy and the RUNJOB data above 10 TeV. (The JACEE helium data are higher.) The 
solid line normalized to the all-particle spectrum at 100 GeV and extrapolated to 1 PeV 
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with an E~ 2 ' 7 power law is consistent with the all-particle measurements of Grigorov 
but about a factor 1.5 above the sum of the RUNJOB spectra at 100 TeV. (Using the 
higher JACEE measurement of helium would bring the emulsion chamber measurements 
into agreement with the E~ 2 ' 7 extrapolation at 100 TeV.) An external motivation for us- 
ing the hard extrapolation from the spectrometer measurements comes from the Super-K 
measurement of the flux of atmospheric neutrinos, where the best fit to the data requires 
an even harder extrapolation of the primary spectrum above 100 GeV 



10]. 



Even with the hard E 2 7 spectrum, the extrapolation of the direct measurements 
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comes in somewhat below the air shower measurements. Moreover, the all-particle spec- 
trum of RUNJOB ^| is somewhat below the lowest-energy air shower data, so there may 
be some systematic offset between the direct measurements and the measurements of air 
showers. There is also a slight correlation between hardness of the fitted spectrum and 
primary mass in this energy range There is, however, no sign of the sharp steep- 
ening of the proton spectrum, below 100 TeV, which would be the case if the maximum 
rigidity accessible in most galactic cosmic-ray accelerators were around 100 TV or below as 
originally estimated for diffusive shock acceleration by expanding supernova remnants 



III. MODELLING THE KNEE 



Although there are systematic differences among the various air shower measurements 



ow 
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in the knee region, all show that the spectrum steepens from —2.7 or slightly harder be. 
1 PeV to —3 above 10 PeV. Figured shows several results. Two of the measurements 
Q] show evidence of structure in the knee region, while the others 0, [zj] do not. 

The figure also shows the spectra of two of the measurements (CASA-MIA and Tibet) 
with the energy assignments of the latter shifted down by a factor of 0.8, which may be 
taken as an indication of systematic uncertainties in energy assignment. The shapes of 
these two spectra agree well with each other. 

Some authors [23j have suggested that a single source contributes significantly to the 
flux of cosmic-rays in the knee region. The possible offset between the air-shower measure- 
ments and the extrapolation of the direct measurements leaves room for such a possibility. 
Others 24,3 

argue, however, that the overall smoothness of the spectrum indicates a sin- 
gle population, with a secondary acceleration mechanism boosting some of the cosmic rays 
accelerated by galactic supernova explosions to higher energy. Several authors j^ . I27I [23] 
approach the problem from the high-energy end by modelling the highest energy spectrum 
as arising from a cosmological distribution of sources and subtracting the extra-galactic 
population from a model of the supernova-accelerated population to see if an extra, high- 
energy galactic source is required to fill in the energy region above the knee ("population 



B" [26]). This is an open question at present. (For a recent review and further references 
see j29l|.) 

In any case it is interesting to estimate the power required to produce the high-energy 
end of the galactic cosmic-ray spectrum. To do so, I assume a model of galactic propagation 
with diffusion characterized by an equivalent leaky box model with a escape time given by 

r esc = 2 x 10 7 yrs x £-°- 33 (2) 

which applies for all energies. With this propagation model (motivated by lack of 
anisotropy at high enery), the source spectrum at low energy must be E~ 2/i7 to give 
the observed E~ 2 ' 7 spectrum inside the galaxy. Assuming a maximum rigidity of 1 PV 
for this low-energy component, with a composition as measured by an emulsion chamber 
experiment ||J, the assumed low-energy component can be subtracted from the observed 
spectrum. The total power required to account for the observed spectrum up to 1 EeV 
is then ~ 2 x 10 39 erg/sec. Such an estimate is obviousy very model dependent, but not 
unreasonably large compared to what might be available in individual galactic sources. 



IV. PRIMARY COMPOSITION FROM AIR SHOWER MEASUREMENTS 

What is needed to make progress is a precise knowledge of the energy-dependence of the 
major groups of nuclei (p, He, CNO, heavy) through the knee region. If there are several 
important groups of sources with different maximum rigidities then there should be a 
corresponding sequence of Peters cycles, perhaps characterized by different compositions. 
The transition to extra-galactic cosmic-rays would be characterized by a transition from 
heavy nuclei (from the highest galactic source) to the light component of a cosmological 
distribution of sources. A recent summary of composition measurements with air showers 
in the knee region is given in the rapporteur paper of Matthews at the last ICRC 3l| (See 
also Swordy et al. |3lj|.) Overall, the evidence suggests a change in composition toward 
hearvier primares at higher energy as expected. However, results of different measurements 
disagree in detail and the picture is not very clear. As noted in the introduction, the 
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best evidence for the sequence of increasingly heavier groups of nuclei comes from the 
KASCADE experiment j|, which uses the ratio of muons to electrons as a probe of 
primary composition. 

Ratio of muons to electrons in the shower front is a classic probe of primary composi- 
tion. At each interaction of the nucleons in an air shower, roughly 1/3 of the energy not 
retained by the projectile nucleon is transferred to the electromagnetic component of the 
cascade via tt° — > 77 and 2/3 to charged pions. The charged pions either reinteract or 
decay depending on their energy and the depth (and hence the density) in the atmosphere 
where they are produced. Charged pions that interact contribute further to the electro- 
magnetic component while those that decay feed the muon component of the air shower. 
Comparing nuclei of the same total energy, charged pions reach the energy at which they 
can decay earlier in the cascade for heavy primaries than for protons because the initial 
energy per nucleon is lower by E tota i/A. As a consequence, the ratio of the muonic to the 
electromagnetic component of an air shower is larger for heavy primaries. 

To measure this ratio with greatest sensitivity requires an observation near shower 
maximum so that the size of the shower is well correlated with the primary energy. In 
general, showers in the knee region and somewhat above are observed after shower maxi- 
mum so that fluctuations in the relation between observed shower size and primary energy 
are important. KASCADE is a surface array near sea level so that these fluctuations are 
large. 

Muons detected in KASCADE typically have energies of a few GeV. An alternative, 
realized with EASTOP-MACRO Q and with SPASE- AMANDA Q, is to sample the 
muon component with a deep underground detector and the electromagnetic component 
with an array on the surface above. Such a setup selects muons with sufficiently high 
energy at production to penetrate to the deep detector. The high-energy muons generally 
come more from the fragmentation region of phase space of the hadronic interactions 
in the shower, whereas than the low-energy muons reflect more the less well-understood 
central region. On the other hand, the multiplicity of high-energy muons is small so that 
fluctuations in muon number are more important. Some of these points are illustrated 
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in Fig. |31 which compares Nu vs. N e for high-energy (> 0.5 TeV) muons and for low- 
energy muons in 10 15 and 10 16 eV showers simulated with CORSIKA. In general, the 
mass resolution is somewhat better for the low energy muons provided that the sampling 
is good enough to get a good measure of the muon number in individual showers, while 
the energy resolution is better for coincident events with high energy muons, particularly 
if the surface array is at high altitude. This complementarity was pointed out by Ralph 
Engel [3J. 
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FIG. 3: Correlation between muons and shower size at the surface. Left panel, muons at 2 km in 
the ice; right panel, muons at the surface. 



It is interesting that analyses of both EASTOP-MACRO [3J and SPASE- 
AMANDA j^] suggest an increase in the mean primary mass in the decade of energy 
above a PeV, that is, through the knee region. As a muon detector, MACRO has the 
advantage of resolving muon tracks, whereas AMANDA reconstructs the light pool gen- 
erated by the bundle, which gives a measure of the energy deposition of the muons. A 
limitation of both experiments is the relatively small sampling area of the deep muon 
detectors, which has to be accounted for in the comparison between data and simulations. 

IceCube, with its surface component IceTop, now under construction at the South 
Pole .35] will have a much larger acceptance. Engineering data from the first year of 
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operation of IceCube with four IceTop stations and one string demonstrate the ability 
of IceCube to reconstruct events with few nanosecond accuracy over distances of more 



than two kilometers 



361 ] . Since January 2006 IceCube has been operating with 16 surface 



stations and 9 strings of detectors in the ice. Each string contains 60 digital optical modules 
(DOMs) evenly spaced in the clear ice between 1450 and 2450 meters below the surface. 
Each IceTop station consists of two ice Cherenkov tanks separate from each other by 10 
meters and located 25 m from the top of the corresponding IceCube string. Each tank is 
instrumented with two of the same DOMs used in the ice. With a nominal grid spacing of 
125 m, the acceptance of the current partial array is ~ 1500 m 2 sr, which allows detection 
of coincident events approaching 10 17 eV. The acceptance of the full IceCube as a three- 
dimensional air shower array will be approximately 0.3 km 2 sr, allowing measurement of 
coincident events up to an EeV. It is scheduled for completion in 2011 and will be operating 
in the meantime as new detectors are added. 



V. TRANSITION TO EXTRA-GALACTIC COSMIC RAYS 

Analysis of measurements of energy-dependence of the depth of shower maximum with 
HiRes |37l ] suggest a change in composition from heavy to light as energy increases from 
10 17 to 10 18 eV. Shower maximum is deeper in the atmosphere for protons than for heavy 
nuclei of the same total energy. The measured depth of shower maximum increases more 
rapidly than expected from model calculations, indicating an increasing fraction of protons 
as energy increases. A similar change was observed in analysis of the Fly's Eye data 38] 
but at a higher energy (above 10 18 eV). It was recognized at the time that such a change 
of composition could be a signature of the transition from a population of galactic cosmic 
rays to an extra-galactic population. In view of the more recent data, the signature of a 
transition appears to occur at a lower energy (below 10 18 eV). A discussion and references 
to cosmological scenarios in which the transition to extra-galactic cosmic rays occurs at 
relatively low energy is given by Hillas j^ . 

Several new experiments are planned or in operation that can explore the energy re- 
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gion from the knee to the EeV region to overlap with the threshold region of the giant 
arrays Auger jsjj and (formerly) AGASA Q. The most advanced of these if KASCADE- 
Grande 4l|, a large sea-level array that includes the original KASCADE as a subarray. 
IceCube is under construction as described in E1IVI Telescope Array [42j and its low-energy 
extension, TALE are under develo pm ent in Utah. There is a propsed array of atmo- 
spheric Cherenkov detectors, TUNKA [4J|, in Russia, and use of the radio technique for 
a large acceptance air shower detector is being explored j^jj. Understanding the transi- 
tion from a galactic to an extra-galactic population of cosmic rays is an interesting and 
important goal for the near future. 
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